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Abstract 

The evolution of single-particle energies with varying isospin asymmetry in the shell 
model is an important issue when predicting changes in the shell structure for exotic 
nuclei. In many cases pseudospin partner levels, that are almost degenerate in energy 
for stable nuclei, are relevant in extracting the size of the shell gaps. A breaking of 
the pseudospin symmetry can affect the size of these gaps and change the magic 
numbers accordingly. The strength of the pseudospin splitting is expected to depend 
in particular on isovector-dependent and tensor contributions to the effective nuclear 
interaction. A description employing supersymmetric quantum mechanics allows to 
derive a pseudospin symmetry breaking potential that is regular in contrast to the 
pseudospin-orbit potential in the conventional relativistic treatment. The derived 
perturbation potential provides a measure to quantify the symmetry breaking and it 
can be employed to improve mean-field calculations in order to better reproduce the 
experimentally observed shell evolution. General potentials with exact pseudospin 
symmetry are obtained that can be used in relativistic mean-field Hamiltonians. 
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1 Introduction 



The single-particle shell model with strong central and spin-orbit potentials 
is a cornerstone of nuclear structure physics since about sixty years. It was 
introduced by Haxel, Jensen and Suefi [T][2][5P] and Goeppert-Mayer |0 in 
order to explain the experimentally observed occurrence of magic numbers for 
nuclei close to the valley of /3-stability. In this picture, nucleons move as inde- 
pendent particles in a spherical mean-field. The central potential usually has 
a shape similar to the nuclear density distribution and the spin-orbit potential 
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is peaked at the nuclear surface. The large spin-orbit potential explains the 
observed energy splittings between two spin-orbit partners with total angular 
momentum j> = / ± 1/2 and parity vr = (—1)' for a given orbital angular 
momentum /, see the left part of figure [H This essential feature of the model 
generates the particular shell structure with the appearance of major shell 
closures at nucleon numbers 2, 8, 20, 28, 50, 82, 126, . . . and subshell clo- 
sures at, e.g., 14, 40, 114. The single-particle wave functions of the mean-field 
calculation can serve as a basis in order to construct fully antisymmetrized 
many-body wave functions. Using a suitably chosen subset in the many-body 
space with an appropriate residual interaction, the excitation spectrum of a 
nucleus can be described with high accuracy by diagonalizing the correspond- 
ing many-body Hamiltonian, see, e.g., ref. [B]. 

Experiments with exotic nuclei suggest a change of the well-known shell struc- 
ture outside the valley of stability and the appearance of new (sub-)shell clo- 
sures often attributed to changes in the spin-orbit splitting or halo effects 
[71l8ll9][T0l[TT][TMMl[TM6^^ It is very important 

to understand the evolution of the single-particle energies with a change in 
isospin asymmetry and to identify the origin of this effect. A change in the 
size of shell gaps can have important consequences in r-process nucleosyn- 
thesis P71l28f29f30ll31] . In recent years, contributions to the nucleon- nucleon 
potential that change the usual shell structure when the isospin asymme- 
try of the system increases have been discussed extensively in the literature 
|32|33|34|35|36|37ll38|39|40|41|42|43ll44j . There are mainly two contributions to 
the nucleon-nucleon potential that have been considered: 1. a central interac- 
tion of Majorana-type VM{r)ai ■ 0^2^^ ■ ^2 (spin-isospin interaction) and 2. a 
non-central isospin-tensor interaction V7T('")>5'i2'ri ■ T2 with the tensor operator 
S12 = 3(cri ■ ri){a2 ■ r2)/{rir2) — ffi ■ 0^2, where the absolute strength of the 
potentials depends on the particular radial functions VM(r) and ViT{r) with 
r = |ri — r2|. The monopole part of these interactions is found to be strongly 
attractive between protons and neutrons when the spins of the two nucleons 
are anti-aligned, i.e., for nlj^ and ulj, (or nlj^ and ulj, ) with j< = I — 1/2, 
j'-^ = I' + 1/2 (in case of the central spin-isospin interaction one has / = /'). 
Hence, the filling of one level with protons (neutrons) will lead to a strong shift 
of the corresponding neutron (proton) states changing the spin-orbit sphttings 
and, depending on the involved levels, the size of the shell gaps. 

Of particular importance is the behaviour of the levels that define the con- 
ventional shell gaps. It is interesting to note that for the (sub-) shell closures 
at 14, 20, 28, 50, 82, 114 and 126 at least one level defining the gap belongs a 
pair of so-called pseudospin partners. These states are marked by solid circles 
in the central level scheme in figure [H The total and orbital angular quantum 
numbers of the partner states are related to the so-called pseudo orbital an- 
gular momentum / by j< = / — 1/2, /< = / — 1 and j> = / + 1/2, /> = / -|- 1, 
respectively. Both levels have the same parity vr = (—1)'+^. In addition, the 
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principal quantum number n of the state with j< = / — 1/2 is larger than the 
principal quantum number of the state with j> = / + 1/2 by one unit. For 
example, the doublets {2si/2, lc?3/2}so = {lpi/25 lp3/2}ps with n = 1, / = 1 and 
{2p3/2, l/5/2}so = {lc?3/2, 1(^5/2 }ps with n = 1, / = 2 are pseudospin partners. 
The indices so and ps denote the spin-orbit and pseudospin representation of 
the states, respectively. Many partner levels of pseudospin pairs are experi- 
mentally found to be almost degenerate in energy suggesting the existence of 
a particular symmetry in the mean-field Hamiltonian for the single-particle 
states. This feature has been noticed already more than forty years ago and 
the concept of pseudospin and pseudo-orbital angular momentum was intro- 
duced [l5P6] . A breaking of the pseudospin symmetry in stable and exotic 
nuclei is potentially important for a change of the level structure and the ap- 
pearance of new magic numbers. It is important to identify the contributions 
to the nuclear interaction that are responsible for a breaking or the restoration 
of the pseudospin symmetry. To this end, a quantitative measure of the pseu- 
dospin splitting is needed that can be used to compare and improve theoretical 
models. 

Bohr, Hamamoto and Mottelson [17] discussed a transformation of the form 
a ■ rjr that connects the angular momentum part of the pseudospin partner 
levels (see eq. ([7])). However, it cannot account for the change of the principal 
quantum number n. Hence, a dynamical symmetry involving the momentum 
operator p is required. Also it was observed that the SU(3) and pseudo-SU(3) 
states of the shell model are related by transformations that are elements of 
the orthosymplectic superalgebra Osp(l/2) [I8]|^ . The many-particle pseudo- 
SU(3) theory is successfully used in the description of heavy deformed nuclei 
in the context of the interacting boson model (IBM) applying supersymme- 
try, see, e.g., ref. [50]. Eventually, the relevant momentum helicity operator 
a ■ p/p was identified that accomplishes the pseudospin transformation and 
it was noticed that the Dirac equation of free or in-medium nucleons (in the 
Dirac-Brueckner approach) is invariant under this transformation [5T][52] . The 
essential step in explaining the pseudospin symmetry was the application of a 
relativistic description. In this approach, e.g., in the form of relativistic mean- 
field (RMF) models, the single-particle states are described by Dirac spinor 
wave functions with large upper and small lower components. Ginocchio ob- 
served that the second-order differential equation for the lower component is 
approximately invariant with respect to the pseudospin symmetry transfor- 
mation in RMF models with large scalar and vector potentials of nearly equal 
magnitude. In addition, the radial wave functions of this lower component are 
almost identical for two pseudospin partner states |53|l54f55f56f57f58f59f60j . 
The pseudospin symmetry was studied in various models for nuclear structure 
without and with deformation |61f62ll63f64f65ll66f67f68ll69ll70f71ll72ll73f74j and 
in nucleon-scattering [75]|76ll77f 781179] . A symmetry breaking pseudospin-orbit 
potential can be derived from the wave equation of the lower component. In 
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principle, this potential should be a measure for the strength of the symmetry 
breaking. It should appear as a small perturbation in comparison with the 
dominating central potential in a mean-field description of nuclei. However, 
this potential exhibits a pole at finite radii and its position depends strongly 
on the level energy itself. Therefore, one cannot really estimate the strength 
of the splitting and test how the choice of the various contributions to the 
interaction influences the symmetry restoration or breaking as observed ex- 
perimentally. The problems related to this peculiarity have been discussed in 
the literature, see, e.g., ref. 

Besides the single-particle levels that constitute pseudospin pairs there are 
several states that are (pseudo-)unpaired in a nucleus. There are all the 
{npi/2}so = {{n — l)si/2}ps states corresponding to / = irrespective of the 
principal quantum number n and the states for given pseudospin I > with 
(pseudo) principal quantum number n = 1 (n = 0). E.g., for the pseudoan- 
gular momentum / = 1 the lsi/2 level has no pseudospin partner whereas for 
n = n > 1 the pairs {{n + l)si/2, n(i3/2}so = {npi/2,nP3/2}ps appear. See the 
right part of figured! The latter feature is a typical observation for a system of 
levels respecting exact supersymmetry. It is therefore tempting to apply super- 
symmetric quantum mechanics to the nuclear single-particle shell model and 
to go beyond the harmonic oscillator case [^8|H9] . As will be shown below, the 
approach employing supersymmetry allows to extract a pseudospin breaking 
potential that is regular in contrast to the standard relativistic description. 
From the condition of vanishing perturbation potential, general potentials can 
be derived that exhibit exact pseudospin symmetry. In order to compare the 
standard description of pseudospin symmetry with the supersymmetric ap- 
proach, the relativistic Hamiltonian will be used as a starting point in the 
theoretical considerations below. In principle, however, the supersymmetric 
approach can also be applied to non-relativistic Hamiltonians. 

The paper is organized as follows. In section [2] the basic equations for describ- 
ing single-particle states in relativistic models are given and the notation for 
quantum numbers, wave functions, potentials and Hamiltonians is established. 
The conventional description of pseudospin symmetry in the relativistic ap- 
proach is discussed in section [3] and a simple example is introduced that will 
accompany the following theoretical presentation. In section HI first the con- 
cept of supersymmetric quantum mechanics is presented as far as necessary to 
solve the problem considered here. More detailed accounts can be found, e.g., 
in refs. [80|81|82|[83] . Then the supermomenta and supersymmetric partner 
Hamiltonians are defined. The form of the reduced supermomenta at small 
radii leads to a particular choice of energy shifts that relate the pseudospin 
partner Hamiltonians of the relativistic upper component wave functions to 
the partner Hamiltonians in the supersymmetric approach. Then, the con- 
nection between the wave functions in the standard relativistic and the su- 
persymmetric description is given. A regular symmetry breaking potential is 
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derived that quantifies tfie pseudospin splitting. It can be used in tfie future to 
compare and improve nuclear structure calculations. In the final subsection, 
general potentials in relativistic Hamiltonians with exact pseudospin symme- 
try are discussed. The last section contains conclusions and an outlook. 



2 Single-particle states in relativistic mean-field models 

The motion of the nucleons in RMF models is governed by strong scalar and 
vector potentials (or self-energies) originating from the minimal coupling of 
scalar and vector mesons to the nucleons. In addition, a tensor potential can 
appear if vector mesons are coupled non-minimally to the meson field tensors, 
i.e. by a derivative coupling. In general, there are both isoscalar and isovec- 
tor contributions to the nucleon self-energies. The details of the particular 
RMF model are not important in the present discussion. See, e.g., references 



to consider only the calculation of the single-particle states for given scalar, 
vector and tensor potentials in the Dirac equation. 

In the following spherical symmetry will be assumed. In this case, single- 
particle levels are characterized by the principal quantum number n = 1, 2, . . ., 
the orbital angular momentum / = 0, 1, . . ., the parity n = (—1)', the total 
angular momentum j = l± 1/2, and its projection Q = —j, + — 
In the relativistic description it is convenient to combine the quantum numbers 
/ and j by introducing the quantity k = ±1, ±2, ±3, . . . where j{K,) = |k| — 1/2 
and 



See table [T]for the relation between the various quantum numbers. Levels with 
different Q for given n and k are degenerate in energy. 

A single-particle state ipnKU with energy Enn (including the rest mass m) of a 
nucleon is obtained by solving the time-independent Dirac equation 




for a more extensive presentation. Here, it suffices 





(2) 



where 



n = a-p + V{r) +(3[m- S{r)] - ipa ■ f{r) 



(3) 
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is the relativistic Hamiltonian with momentum operator p and Dirac matrices 
d and j3 in standard notation. The system of units is chosen such that h = 
c = 1. The scalar, vector and tensor potentials are denoted by S{r), V{r) and 
T(r) = T{r)r/r. Here, the sign convention with a positive scalar field S is 
used. For spherical symmetry it is customary to write the Dirac spinor as 

1 ( FUr)ynn{r) . 
^n,^n = - (4) 



r 



with radial wave functions Fnnir) = Fni^ (r) and GnK.{r) = Gni^ (r) in the upper 
and lower components, respectively, that are normalized according to 



dr 



i^„«(r)|' + |G„.(r)|' =1. (5) 







The spin and angular dependence is contained in the spinor spherical harmon- 



ics 



mims 



where the orbital angular momentum /(/t) of the spherical harmonic Yjm; is 
coupled with the nucleon spin s = 1/2 of the spinor Xsms to the total angular 
momentum j{K). The upper and lower component in (j3]) have an opposite sign 
for the quantum number k, in the spinor spherical harmonics and consequently 
opposite parity with a difference of one in the orbital angular momentum. By 
applying the Hamiltonian ([3]) to the Dirac spinor (jl]) and using the relations 

a--y.n{r) = -y-.n{r) (7) 
r 

and 

a ■ p^y.nif) = -(^--] f(^)y-^n(^) ' (8) 
r r \ar r J 

the set of coupled first-order differential equations 



{Er,,-V + m- S) Gn.{r) =(^-^-'l-T^ F„«(r) (9) 
(K« - 1/ - m + S) FUr) = - l^i- + ^ + T j G„«(r) (10) 
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for the radial wave functions is obtained. A simple decoupling leads to the 
independent wave equations 



HF{K)Fn 
HG{K)Gn 



E F 



with Schrodinger-like Hamiltonians 



:i2) 



_ 1 



and 



' d K ^\ 1 f d K 



T 



V -S + m 



d"^ k,(k, 
+ ^— 



dr"^ 

+V - S + m 



+ 2T- + + — 

r A \ ar r 



T 



(13) 



1 

A 



d^ k{k+1) B' f d K ^ 



dr'^ r 
+V - S + Tn, 



B \dr 



(14) 



respectively, that explicitly depend on the quantum number k. A prime de- 
notes the derivative with respect to the radius r. The quantities 



A{r) = En. - V{r) + m - S{r) 



(15) 



and 



B{r) = En.-V{r)-m + S{r) 



(16) 



depend on the single-particles energy Enn = Enij appearing as the eigenvalue 
in equations (fTTj) and (fT2l) . The Hamiltonians Hp{k) and Hg{k) have the 
same spectrum, even though the potentials have a different form. A{r) can 
be considered as twice the position depending effective nucleon mass rries. It 
is always positive in standard RMF models. The dependence on the energy 
EnK. is a relativistic effect. Since En. contains the rest mass m of the nucleon, 
the variation of A{r) with a change in the nucleon binding energy is very 
small. For vanishing scalar and vector fields S{r) and V{r) at large radii r, 
B{r) reduces to the non- relativistic energy = — m < of the bound 
nucleon. The Hamiltonian (fTTI) for the upper component wave function Fnnif) 
can be written in the familiar non-relativistic form 



d 



+ ^-^^ + K(r) + Vso{r){l- + m 



dr 2meff dr 2meff 



(17) 
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with central and spin-orbit potentials 



VAr) = V-S 



Vso{r) 



V -s 



T 




1 








+ — 




I 




Ar 






2T + 



d 



2mefj dr 



{V + S) 



(19) 



since n{n — 1) = /(/ + 1) and k — 1 = 2(/ ■ s). The spin-orbit potential contains 
contributions from the tensor potential T and the derivative of the sum 1^ + 5* 
whereas the main contribution to the central potential arises from the differ- 
ence V — S. With the typical strong vector and scalar potentials of almost 
similar magnitude inside nuclei, the usual depth of the central potential and 
the strong surface-peaked spin-orbit potential are nicely explained in the rel- 
ativistic approach. A tensor potential is usually omitted in traditional RMF 
calculations but it has a distinct effect on the spin-orbit splitting that can be 
relevant for the size of the shell gaps. 



3 Relativistic description of pseudospin symmetry 



The basic idea in the relativistic description of pseudospin symmetry is a 
comparison of the Hamiltonians Hg{k), Hg{k,') and of the radial wave func- 
tions G(„+i)K(r), GriK'ii^) for the two pseudospin partner levels of the doublet 
{(n -|- l)K,nK'}so- It is easy to see that for every state with k > the pseu- 
dospin partner state has the quantum number k' = —n — l. The corresponding 
orbital angular momenta are I = k — 1 and /' = — (— k — 1) = k, + 1, respec- 
tively. Similarly, the total angular momenta are given by j = k — 1/2 and 
j' = I — K — 1| — 1/2 = K -|- 1/2, respectively. Hence, the relation between / and 
K is given by 

K if K > 



(20) 
-K - 1 if K < 



The single-particle states in the relativistic model with k = —1, i.e. the pi/2 
states, have no pseudospin partner. In this case / = 0. See table [1] for a 
summary of the relations for the quantum numbers. 

A simple example shows the main features in the relativistic description of 
pseudospin partner states. Here, the lowest pair of levels with / = 1, i.e.. 



8 



{2si/2, ld3/2}so, will be considered. The shape of scalar and vector potentials in 
RMF models is often well approximated by a Woods-Saxon form, i.e., S{r) = 
Sofir,R,a) and V{r) = Vofir,R,a) with f{r,R,a) = 1 / {1 + exp[{r - R) / a]} . 
Assuming a radius of i? = 3.8 fm, a diffuseness of a = 0.65 fm and absolute 
magnitudes of 5*0 = 450 MeV and Vq = 370 MeV, respectively, (corresponding 
to the the mean-field potentials of neutrons in a nucleus like ^'^Ca) one finds 
single-particles energies of £^251/2 — m = —15.604 MeV and -^1^3/2 — m = 
-15.424 MeV with the Idg/a level lying only 0.180 MeV above the 2si/2 level. 
The corresponding wave functions of the upper and of the lower component 
in the Dirac spinor are depicted in the top and bottom panel of figure [2J 
The wave functions F2s^^2(^) -^1^3/2 very different where the former 

has an additional node because there exists a stronger bound lsi/2 state at 
-^1^1/2 ~ ''^ = —56.848 MeV without node (except at r = 0). In contrast, 
the wave functions G2si/2(^) and Gid^^^{r) resemble each other closely, a sign 
for the almost perfect pseudospin degeneracy of the two states. Due to the 
different orbital angular momenta and corresponding centrifugal potentials in 
Hf{k) and Hf{h') there is a clear difference in the radial dependence of the 
upper component wave functions at small r. But for the lower component 
radial wave function the r dependence is very similar. 

In fact, the orbital angular momenta in the Hamiltonians Hg{h) and Hg{—h — 
1) of a pseudospin pair [(n -|- 1)k, n(— k — l)]so for k = / > are identical since 
= k = 1{k + 1). The centrifugal potential k(k + l)/(Ar^) in f[T^ is 
invariant with respect to the replacement k — > —n— 1. The difference between 
the potentials 

A([, r) = Ha{K) - Ha{-K - 1) = + 2r^ (21) 



could be considered as a measure for the pseudospin splitting because for 
A(/) = the Hamiltonians Hg{h) and Hg{—h— 1) are identical and yield the 
same energy eigenvalues. 

However, there is a problem in this approach since the quantity B{r) has a 
zero at finite r inside the nucleus with usual scalar and vector potentials in 
RMF models. This is most easily seen in figure [3] where the potential difference 
V{r) — S{r) and the non-relativistic energy -E2S1/2 ~ ^ plotted. They are 
identical at a certain radius r and the quantity f|T6|l is zero. Consequently, 
the difference potential fl2Tl) . shown as a red solid line in figure [3|, has a pole. 
Its position is strongly energy dependent. Thus, A(Z) cannot be considered 
small as compared to the potential V — S in the differential equation and it 
is therefore not a good measure for the symmetry breaking. 

For vanishing tensor potential T = 0, the pseudosp in-orbit potential (pT!) is 
zero for 5' = or S{r) = V{r) + C with a constant C. As a consequence of 
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this condition, the main contribution — S* to the central potential (fT8|) is 
just an overall shift of the energy scale for exact pseudospin symmetry. But 
then there will be no sufficiently strong central potential to bind the states as 
required in a nucleus. 

A case where exact pseudospin symmetry is found is the relativistic harmonic 
oscillator with V{r) = —S{r) = muj'^r'^/A and T = [711l92|93|94ll95] . The 
effective mass meff = A/2 = {E + m)/2 is a constant, B{r) = E — m — mu'^r^ /2 
but B'{r) 7^ and A(/,r) ^ 0. Hence, the identification of A(/,r) as the 
relevant quantity to measure the pseudospin splitting is doubtful. It would be 
convenient to have an independent measure of the symmetry breaking that 
has a regular behaviour for all radii and that vanishes for exact pseudospin 
symmetry. 



4 Supersymmetric description of pseudospin symmetry 

The problems related to the conventional relativistic explanation of pseudospin 
symmetry call for an alternative approach that allows to derive a regular 
symmetry breaking potential which is a small perturbation without a pole 
and useful for a quantitative comparison of different models. The main idea of 
the supersymmetric approach to describe the pseudospin symmetry in nuclei 
is similar to the standard relativistic description. In both cases one does not 
compare the original model Hamiltonians of the pseudospin partner systems 
but closely related Hamilton operators with the same spectra. In fact, it will be 
shown below that a comparison of the supersymmetric partner Hamiltonians 
provides a way to define a proper pseudospin splitting potential. These partner 
potentials are derived from the Hamiltonians Hf{k) and Hp[k') for the upper 
component wave functions of the pseudospin partners k and k' = —k — 1, 
respectively. To this end, the Hamiltonian (fT3l) is written as 

Hf{k) = Hi{k) + E{k) (22) 

with a Hamiltonian Hi{k) of a supersymmetric pair and an energy shift E{k,) 
to be determined below. A corresponding relation holds for the pseudospin 
partner Hamiltonian Hf{k'). 

4-1 General features of supersymmetric quantum mechanics 

It is well known, that every second-order Hamiltonian can be factorized in 
a product of two Hermitian conjugate first-order operators [55|97II98] . In the 
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present case, operators and = [-B^]^ are introduced such that 

H,{k) = BtB- . (23) 

Then, the Hermitian operators 

( ^ B+\ ( -iB-^ , 

Qi{^)= Q,{^) = tQi{^)T = I (24) 

\B- ) [zB- 

are formed that are so-called supercharges with respect to the involution 

r = rt = ^ ° I = 1 (25) 



because {Qi,t} = {Q2,t} = 0. In the next step, the supersymmetric Hamil- 
tonian 

Hs{k) = mK)f = mn)]' =\ ^ I (26) 

H2{k) 



is obtained with the supersymmetric partner Hamiltonians (]23ll and 

H^in) = B-B+ (27) 

on the diagonal. Obviously, Hs{n) commutes with Qi{k) and Q2{i^), i-e., 

[Hs{^),Qi{^)] = [Hs{k),Q2{k)] = (28) 

and the supercharges Qi{K) {i = 1,2) are generators of the symmetry trans- 
formation. The supercharges ( l24l) and the Hamiltonian (l26ll with the commu- 
tators fl28l) and the anticommutator 



{Qi(/^),Q2(/«)} = (29) 

are the most simple example of a supersymmetric algebra. Since Hs{k) is 
the square of the Hermitian operators Qi{K,), all eigenvalues Es{nn) of the 
eigenvalue equation 

Hs{K)'^s{nK) = Es{nK)'^s{nK) (30) 
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with the two-component wave function 
( , 




(31) 



are non-negative. It is easily seen by applying Qiin) to equation fl5U]) . that 
Hi{k,) and H2{k) have the same spectrum of positive energies Es{nn) > 
where the operators B'^ and B~ connect the components of the wave function 
(13T1) by the transformations 



-.ipiiriK) ipiiriK) = 



Es{nK) 



■■i'2inK) 



(32) 



If there is an eigenstate \E'5'(0/t) (n = 0) with energy Es{Ok) = 0, the su- 
persymmetry is called exact because (5i(«;)^5(0K) = and this ground state 
obeys the symmetry of the Hamiltonian Hs{n)- In this case, 



B-iJi{OK) = ij2{0K) = 



(33) 



or 



B+^lJ2{0K) = ^/'i(Ok) = 



(34) 



i.e., the Hamiltonian Hi{n) or H2{k) has an additional state at zero energy 
that is not appearing for its supersymmetric partner Hamiltonian. In this 
paper, the usual convention is chosen such that Hi{k) has a ground state at 
zero energy if the supersymmetry is exact. If there is no state \I'5(0k) at energy 
£^5 (Ok) = 0, the supersymmetry is called broken since Qi{K)'^ s{nK) 7^ for 
all n = 1, 2, . . . and the partner Hamiltonians Hi{k) and H2{k) have identical 
spectra. 



4.2 The supermomentum and supersymmetric partner Hamiltonians 



The main task in the supersymmetric description of pseudospin symmetry 
is the construction of the operators B'^ and B~ . The particular form of the 
Hamiltonian (fT3l) suggests to use the ansatz 



B- 



AJr) 



d 
dr 



Bt 



d 

dr 



(35) 
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with (the principal quantum number n is suppressed in the following if not 
exphcitly needed) 

A,{r) = E,- V{r) + m- S{r) (36) 



and a quantity Qn{r) that will be called supermomentum by dimensional 
reasons. The traditional superpotential is given by 

W.{r) = ^I^^Q^ir) , (37) 



however, it has the dimension energy^/^. The supersymmetric partner Hamil- 
tonians are found as 



1 



Qi - q: 




dr 



(38) 



and 



Ho(k) 



B-Bl = — 



d 



A' d 



Qk + Qk to + /I j„ + 



A'' 



I \2' 



dr^ A, dr 2A, 



AAl 



(39) 



with the symmetry 

H,{k,Q^) = H,{k,-Q^-^). (40) 



If the supersjTiimetry is exact, the supermomentum can be found from the 
ground state wave function '?/'i(0/t) of the Hamiltonian Hi{k) as 

Qn{r) = -^\niJi{^K) (41) 
dr 



because of equation (l33ll . i.e., it is given by the negative logarithmic derivative 
of the groundstate wave function. Comparing equation fl38l) via (l22l) with 
equation ffTTj) leads to the defining equation for the supermomentum Qnif) in 
general. However, it is useful to split off the dependence on k and the tensor 
potential T first and to introduce the reduced supermomentum 

q.{r) = Q^{r) + - + T{r) (42) 



with the results 
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and 



1 



"^"^ +ql-2q^-q'^-2q.T + T^ (43) 



r 



-2T- + r + -^iq^ + - T 

r \ dr r , 



1 

A. 



k{k + 1) 



+ 



K 



+ ql-2q,- + q',-2q^T + T' 



r 



A' d A" 3 (AM 2 
r~ Xrf^ 2X~ 4^2 



These two Hamiltonians show the symmetry 

A' 



ifi(K, g^, T) = H2{-K, -q^ 



2A^ 



(44) 



(45) 



however, it is irrelevant for the present considerations. Subsequently, the defin- 
ing equation for the reduced supermomentum 



ql- 2- + 2T--^]q^-q: = -A^B^-T' (46) 




with 



B^{r) = E{k) - V{r) + S{r) - m 



(47) 



is obtained. Note that Bj^r) depends on the energy shift E(^k) and ^^(r) on the 
eigenvalue E^,. Equation (H^ is the central relation that connects the reduced 
supermomentum ^^(r) with the relativistic potentials S{r\ Vir) and T(r). 
In general, the reduced supermomentum q^ir) is not uniquely determined by 
equation fH6l) and suitable boundary conditions have to be specified. In the 
present case, it will be required that qjyr^ = for r = 0. 

It is convenient, e.g., in numerical applications, to transform the Riccati equa- 
tion fB6l) . an inhomogeneous nonlinear first-order differential equation, into a 
linear second-order differential equation to determine the reduced supermo- 
mentum q,^{r) for given scalar, vector and tensor potentials S{r\ Vir^ and 
T(r). To this end, a new function yj^r) 7^ is introduced that is related via 

= (48) 
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to the reduced supermomentum qK^r). With 

q' =-i^ + ^ (49) 



one immediately obtains the homogeneous differential equation 

y: + (2^ + 2T-^]y'^+ (A^B^ + T') = (50) 



that can be solved with standard techniques. 

Using equation (H6i) the supersymmetric partner Hamiltonian (jUl) assumes 
the form 



(f ^ , A'^ ( d\ , ^, 



dr"^ ^ Ak \ dr 

B 



r 2A, AAl 



It will be useful for the comparison with the Hamiltonian H2{k') of the pseu- 
dospin partner system. 



4.3 Asymptotics of the supermomenta 



The asymptotic behaviour of the reduced supermomentum for r — and 
r — >■ 00 is easily derived from equation 0461) assuming regular qK{r) and regular 
potentials S{r), V{r) and T(r). The reduced superpotential has to vanish for 
r —>■ otherwise a diverging potential would appear in the Hamiltonian Hf{k,). 
For small radii the reduced supermomentum is given by 

g«(r) = — — r + 0(r ) , (52) 



i.e., it increases almost linearly with the radius. At large radii it becomes a 
constant 



g,(r) = V (^. + m) [m - E{k)] + 0{r-') (53) 

for potentials S{r), V{r) and T(r) that approach zero at large radii r. The 
energy shift E{k) with the condition E{k) < m has still to be determined. 
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The asymptotic form (l53l) translates to a dependence 



y^ir) exp 



-^J + m) [m - -^(/t)] r 



} 



(54) 



of the auxihary function in equation fl50l) for r ^ oo. Starting at large 
radii and integrating inwards a stable solution for 2/k('") is found. Of course, 
if the relativistic potentials do not approach zero at large radii, a different 
asymptotic behaviour of the reduced supermomentum will be obtained. 

Considering the relation P2|) . the full supermomentum Qi^{r) approaches the 
same constant as qK{r) for large r. In contrast, the asymptotic behaviour of 
QK.{r) for r — i> is determined by the angular momentum term n/r. For /t > 0, 
one finds \mir^oQi^{r) = — oo and limr^oQ^'ir) = oo for k' = —k — 1 < 0. 
The standard superpotential (1371) shows the same asymptotic limits as the 
full supermomentum. From general considerations of supersymmetric quan- 
tum mechanics it is known, that the type of supersymmetry is determined by 
the asymptotics of the superpotential. If there is a change of sign in Wi^{r) 
when comparing the limits r — with r — > cxd, exact supersymmetry fol- 
lows and a single non-degenerate state at zero energy exists. Conversely, if the 
superpotential Wf^{r) does not change the sign in the above limits the super- 
symmetry is broken and all eigenstates are doubly degenerate with positive 
energy. 

4-4 Application of supersymmetry to pseudospin levels and energy shifts 

With the foregoing considerations, the supersymmetric description can be ap- 
plied to the system of pseudospin levels appearing in nuclei. For given pseudo 
orbital angular momentum / = there is an unpaired state with total 

angular momentum j = Z — 1/2, orbital angular momentum 1 = 1 — 1 and 
energy Ei^, cf. the left part of figure HI Choosing the energy shift as 



the ground state energy of the Hamiltonian Hi{k) = Hp{k)—E{k) is placed at 
zero energy and the case of exact supersymmetry is obtained. The supersym- 
metric partner Hamiltonian H2{k) has the same energy eigenvalues as Hi{k) 
except for the zero-energy level. The orbital angular momentum in H2{k) is 
I = K = I as seen in the centrifugal potential in equation (jH]). 

The pseudospin partner levels with j' = / + 1/2 and I' = I + 1 are eigenstates 
of the Hamiltonian Hf{k') with k' = —I — 1 < 0, see the right part of figure 
m The lowest state has a larger energy Ei,^! than the ground state of Hp{k). 



E{k) = E-_ 



Ik 



(55) 
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For this system the case of broken supersymmetry apphes where Hf{k,') = 
Hi{k') + E{k') and the supersymmetric partner Hamiltonian H2{k') + E{k') 
have identical eigenvalues. The orbital angular momentum in H2{k,') is again 
given by /' = K = / due to + 1) = (— k — — 1 + 1) = k{k + 1). The 

energy shift E{k') is now chosen such that the asymptotic form of the reduced 
supermomentum qK'ir) is the same as for the reduced supermomentum ^^(r) 
in the limit r — > 0. This condition leads with equation fl52l) to the relation 

^,(0)5,(0) + T'(0) = - [A^iO)B^iO) + T'(0)] (56) 



since 2n' + 1 = —(2k + 1). For exact pseudospin symmetry, E(^n+i)K = Enn' 
and consequently Ai^{r) = A^'ir). Even if the pseudospin symmetry is broken, 
and Af^i are almost identical and the difference can be neglected as a small 
relativistic correction. Without the contribution of the tensor potential T'(0), 
the energy shift E{k,') is determined by 

E{k') -m + S{0) - V{0) = - [E{k) - m + S{0) - V{0)] . (57) 



In the top panel of figure O the reduced supermomenta qK^r) and q^'ir) of 
the pseudospin pair levels with k = 1 and k' = —2 are shown for the ex- 
ample in section [3l With the present choice of the energy shifts E{1) — m = 
—56.849 MeV and E{—2) — m = —102.680 MeV, the reduced supermomenta 
display the same approximately linear behaviour for r ^ as required. In fact, 
they are almost identical for radii below 5 fm. At large r they become constant 
and approach different values for r — > cxd as predicted by equation fl53|) with 
qi{r) < q-2{,f) because E{1) > E{—2). The full supermomenta Qi{r) and 
Q-2{r) are shown in the bottom panel of figure [51 They contain the angular 
momentum contribution and, hence, they diverge for r —>■ 0. The change of 
sign for Qi(r) indicates the case of exact supersymmetry with a single state 
at zero energy for Hi{l). In contrast, there is no change of sign for Q_2(r) 
corresponding to the case of broken supersymmetry. 



4-5 Relation of wave functions 



The supersymmetric description of pseudospin degenerate states allows to 
establish relations between the wave functions in the Dirac spinor (jlj) and 
in the supersymmetric state vector (1311) . According to equation the lower 
component wave function of the spinor is given by 

Gn.{r) = ^{^---T] FUr) (5J 
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for the pseudospin partner state tik with k > and similarly for the state n'n' 
with n' = n — 1 and n' = —k, — 1. Using equation (132!) the supersymmetric 
partner 

V'2(nK, r) = . IqnK - - - T + ) i>i{nK, r) (59) 



of the upper component wave function ipiijin^r) = Fn^ir) is obtained where 
the eigenenergy Esinn) in the supersymmetric formulation is identified with 
the energy E^k — E{n) = En^ — -Eqk- With equation fl58l) one finds 



Er, 



E{k) 



-Ar, 



-FnJr) + Gr, 



(60) 



and 



En'n' - E{k') 



(61) 



for the two pseudospin partner states. These functions are normaUzed as 

oo CxD 

y" dr |?/'2(nK,r)|^ ^ J \FnKir)f (62) 





and correspondingly for the second pair. In the case of exact pseudospin sym- 
metry, the wave functions ip2{nK, r) and ip2{n'K,', r) are identical provided that 
the upper component wave functions Fnnir) and Fn'K'ij) in the Dirac spinors 
are normalized to the same value. This is beautifully seen in figure [6]for the ex- 
ample of section [3] with a much better agreement of the two pseudospin partner 
wave functions than in the standard approach, cf. figure [2J In contrast to the 
wave functions Gnn{r) and (r), there is no zero in the functions ijj2{nK, r) 
and il)2{n'K\ r) because they are the lowest states of the supersymmetric part- 
ner Hamiltonians i?2(^) and H2{k') with the same orbital angular momentum 
/ = [ = K = -k' - 1 = /'. The factors g„K(r)/yl„«(r) and g„'K'(r)/A„/«/(r) in 
equations fIBUl) and flMl) are very small at small radii r explaining the obser- 
vation of the close similarity of the wave functions GnK(r) and Gn'K,'{T) in the 
usual comparison. At larger radii, however, the contribution from the upper 
component wave functions Fnnij") and Fn'K'ij), respectively, becomes more 
important. 
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4-6 Pseudospin symmetry breaking potential 

The Hamiltonians H2{k) + E{k,) and H2{k,') + E{k,') of the supersymmetric 
partners for the pseudospin partner levels are almost identical. The difference 
is given by the potential 

As{l, r) = H^i^) + E{k) - [H^W) + E{k')] (63) 

\fA dr ^/A Ar 

with 1 = h = —k' — 1 and A = A^^ where equation fIFI]) has been used and 
the difference between and has been neglected. The quantity As{l.,r) 
is a regular function for all radii r with the same contribution from the tensor 
potential T(r) as the previously defined pseudospin-orbit potential A(/,r) in 
(l2T]) . However, the pseudospin breaking potential (!63|) is a regular function 
for all r without a pole as observed for A(/,r). Without tensor potential, the 
pseudospin symmetry breaking potential depends on the derivative containing 
the reduced supermomenta of the pseudospin partner systems instead on the 
derivative of the difference V — 5* as in the traditional approach fl2T|) . From 
figure [3] it is evident that As{l.,r) with I = 1 for the example in section 
[3] is really small as compared to V{r) — S{r) and can be considered as a 
perturbation. In a relativistic model with tensor potential T{r), both the non- 
relativistic spin-orbit potential flTI?]) and the pseudospin symmetry breaking 
potential (163!) are directly affected by T(r) with the same functional form, 
only the prefactors n — 1 and 2k + 1 are different. The dependence of As(/, r) 
on the scalar and vector potentials S{r) and V{r) is, however, more subtle 
through the Riccati equation for the reduced supermomenta. 



4-7 Potentials in relativistic Hamiltonians with exact pseudospin symmetry 

In general, there are three potentials, S{r), V{r) and T{r), in the relativistic 
Hamiltonian With one constraint As = to achieve exact pseudospin 
symmetry, only two independent functions of these potentials can be specified 
freely (subject to some conditions). Since the difference potential fl63|) is a 
function of A{r) (containing the sum V{r) + S{r)) and T(r), it is reasonable 
to use these quantities as independent variables and to determine the reduced 
supermomenta, q^ir) and q^'ir) and finally the scalar and vector potentials, 
S{r) and V{r). This approach will be followed in the considerations below. 
As an alternative, the potentials in the relativistic Hamiltonian could be de- 
rived from two given supermomenta qKir) and qn'if)- In order to simplify the 
notation, the quantities 
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UJ ■ 



En 



2k + 1 
E(k) + E(k') 



(64) 
(65) 



are introduced replacing the energy shifts E{k,) and E{k,'). 



From the condition As{l,r) = one obtains the general solution for the 
difference 



q^{r)-q,ir) = -{2K+l)C{r) 
with the function 



(66) 



T(r') dr' 



(67) 



where D is an integration constant. Writing 



(lK{.r) = q{r) + T(r) - ^ C{r) 

2k + 1 

g^,(r) = g(r) + T(r) + — — C(r) 



(68) 
(69) 



the difference of the defining equations fH6l) for qK.{r) and gK'(r) leads to the 
solution 



q{r) 



{Auj - C')r - C 
1 + Cr 



Aur 



d 



1 + Cr dr 



ln(l + Cr) 



(70) 



for the auxiliary momentum q{r). The sum of the two defining equations gives 
in combination with the definition of A{r) the scalar and vector potentials 



S{r) = ^ - - A{r) - U{r)] + m 
V{r) = 1 + - A{r) + U{r)] 
with the auxiliary potential 



(71) 
(72) 



U{r) = V{r)-S{r)+m-Efi 



(73) 



1 



+ T) g-T 



1 A^ 

r^ A 



q' + {2k + 1] 









'9. 
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that is invariant under the replacement 1. For the determi- 

nation of the potentials S{r) and V{r) in the relativistic Hamiltonian that 
exhibits exact pseudospin symmetry, the given functions A{r) and T(r) have 
to be regular, but not every choice is admissible. At least, one has to require 
that A{r) > for r > and that T(r) is selected such that the integral 
in equation fl67j) is finite for all r. In addition, there is the constraint that 
C(r) has to be zero at r = because ^^(O) = qn'iO) = by construction of 
the supersymmetric partner Hamiltonians and that C{r) > 0, otherwise the 
momentum q{r) is not regular for all radii. It is worthwhile to note that the 
potential (1751) depends on the quantum number k if C{r) ^ and that exact 
pseudospin symmetry is found only in the particular doublet of states with 
given K, = I = K, and /t' = — / — 1 if the relativistic potentials are assumed to 
be identical for all states. 

It is instructing to discuss the resulting potentials for a few particular choices 
of A{r) and T(r). There are two major cases depending on the value of A{r) 
at r = 0. 

(1) A{0) = but A{r) > if r > 0. In this case, which might not be realized 
in a physical system, the integration constant D in equation fl671) can 
have any non- negative value as long as C{r) > for all r. If C{r) is not 
identically zero, the reduced supermomenta of the pseudospin partner 
system with different I are different and the potential U{r) explicitly 
depends on k. This can occur even for vanishing tensor potential if D > 0. 

(2) A{0) 7^ 0, the case usually encountered in physical systems. It follows 
that D = and that C{r) ^ only if there is a tensor interaction. See 
the discussion in the paragraph above. 

In most relativistic Hamiltonians without a tensor interaction and everywhere 
finite effective mass of the nucleon there is still a large variation of potentials 
that exhibit exact pseudospin symmetry. For T(r) = and then C(r) = for 
all radii r, the reduced supermomenta of the pseudospin partners are identical 
with 

A(r) 

qnir) = q^c'ir) = q{r) = -i^^r (74) 



and the auxiliary potential becomes 
U{r 



1 



r A[r) 



4 



(75) 



i.e., a modified harmonic oscillator potential with position depending effective 
mass. With the scalar and vector potentials (!7T|) and (!72l) . respectively, the 
central and and spin-orbit potentials (ITSl) and (fT9l) are found as 
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V;(r) = E,-m + U{r) + ^Vso{r) (76) 

in the non-relativistic Hamiltonian ( |T7I) . The standard harmonic oscillator is 
obtained for constant A{r) = 2m, T(r) = and E{k) = m + (2k + 1)07/2 such 
that Eq = m. Then the scalar and vector potentials are 

S{r) = -V{r) = -^c^V (78) 



leading to the non-relativistic central and spin-orbit potentials 

K(r) = -uj^r^ Vsoir) = . (79) 



In general, the choice of the position depending effective mass mefT(r) = 
Af^{r)/2 has an important effect on the potentials that appear in the non- 
relativistic Hamiltonian. However, a reasonable effective mass should approach 
the free nucleon mass m for large radii r. Then the potentials will approach 
the harmonic oscillator form in the limit r —>■ oo. 



5 Conclusions and outlook 



The occurrence of almost degenerate pseudospin partner levels in atomic nuclei 
is usually attributed to a symmetry of the underlying relativistic Hamiltonian 
as it appears, e.g., in RMF models. The lower-component radial wave func- 
tions in the Dirac spinors are found to be very similar. From the differential 
equation of these wave functions a pseudospin-orbit potential is extracted that 
breaks the pseudospin symmetry. However, it is not a regular function and it 
cannot be considered as a small perturbation. In addition, this potential does 
not vanish for the usual harmonic oscillator potential which exhibits exact 
pseudospin symmetry. 

In the supersymmetric approach to pseudospin symmetry, the supersymmetric 
partner Hamiltonians of the upper-component wave functions are compared. 
It leads to a regular symmetry breaking potential that is small in compari- 
son with the dominating central potential in the mean-field Hamiltonian. The 
corresponding partner wave functions are identical in the case of exact pseu- 
dospin symmetry. Employing both exact and broken supersymmetry, all levels 
in a sequence of states for given pseudo orbital angular momentum I > 0, in- 
cluding the lowest unpaired state, fit into the scheme. This lowest level defines 
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the reduced supermomentum and consequently the potentials in the Hamil- 
tonians for both pseudospin partner systems in the case of exact pseudospin 
symmetry. All higher lying pairs are degenerate in energy. From the condition 
of vanishing symmetry breaking potential, general relativistic potentials that 
exhibit exact pseudospin symmetry can be derived. Modified harmonic oscil- 
lator potentials with position depending effective mass are the most simple 
class. 

The present description relies on the differential equations for the radial wave 
functions and the corresponding first-order differential operators in the factor- 
ization of the Hamiltonian for given quantum number k. It is also possible to 
generalize the approach to a form independent of a particular k. In this case 
the operators 



a 



T 



1 



A 



B- 



1 



a - (f + 



in the generators ( 12^ of the symmetry transformation act on the full wave 
functions in coordinate space. These and more formal aspects of the super- 
symmetric approach can be studied in the future. 

The symmetry breaking potential depends on derivatives of the supermomenta 
for the two pseudospin partner states and the tensor potential. The supermo- 
menta itself are determined via a Riccati equation by the scalar, vector and 
tensor potentials in the relativistic Hamiltonian. A simple example shows that 
the reduced supermomenta of pseudospin partner states are very similar rising 
approximately linear at small radii as in the case of the harmonic oscillator. It 
seems that all saturating potentials with a fiat bottom will show an approxi- 
mate pseudospin symmetry. This question can be explored in more detail by 
expanding the reduced supermomenta in a power series in terms of the radius 
r. 

Calculating the supermomenta for actual selfconsistent RMF models will help 
to quantify the various sizes of the observed pseudospin symmetry breaking. 
The study of the symmetry breaking potential will eventually lead to improved 
RMF parametrizations since it can help to identify the relevant modifications 
of the model that are required to obtain a better description of the exper- 
imentally observed breaking and restoration of the pseudospin symmetry in 
nuclei. It will be interesting to see how the isospin dependence of the effective 
interaction changes the degree of pseudospin symmetry breaking. Of special 
importance is the tensor interaction (neglected in most RMF models) that 
affects both the spin-orbit and the pseudospin-orbit splitting. 

The method that is developed in this paper is not restricted to a relativistic 
description of the single-particle states. In principle, it can also be applied 
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to the usual Schrodinger equation for the single-particle wave functions as it 
appears, e.g., in the Skyrmc Hartrcc-Fock method, because it relies only on 
the Schrodinger-equivalent differential equation of the upper component in the 
Dirac spinor. Instead of the scalar, vector and tensor potentials in the rela- 
tivistic Hamiltonian, there are also at least three independent quantities that 
characterize the mean-field in the non-relativistic Hamiltonian: the position 
depending effective mass, the central and the spin-orbit potential. The relation 
between these functions and the reduced supermomentum will be different as 
compared to the relativistic description, however, the same principles of the 
supersymmetric method can be applied. 
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Note added in proof 

Pseudospin symmetry was found as a particular limit in the application of 
supersymmetry to the Dirac equation in Refs. [Lev04,Lev05] by directly fac- 
torizing the Dirac Hamiltonian. 
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Tables 



Table 1 

Relation between the various quantum numbers in the relativistic description of 
single-particle states. Pseudospin partners with the same pseudoangular momentum 
I are placed in the same column. 
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Fig. 1. Single-particles states in the conventional shell model description with 
spin-orbit partners j = / it 1/2 for given orbital angular momentum I (left) and 
in the supersymmetric approach with pseudospin partners j = / it 1/2 for given 
pseudo-angular momentum I (right). The pseudospin partner states are marked 
with a solid circle in the central level scheme. Additionally the magic numbers for 
shell closures are given. Subshell closures are enclosed in parentheses. 
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Fig. 2. Wave functions of the upper (top panel) and lower (bottom panel) compo- 
nent in the Dirac spinor for the 25^/2 (^^^d solid line) and 1^3/2 (blue dashed line) 
pseudospin partner levels. 
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Fig. 3. Pseudospin symmetry breaking potentials in the standard relativistic de- 
scription (red solid line) and the supersymmetric approach (blue dashed line). The 
potential difference V{r) — S{r) is denoted by a black dotted line and the energy of 
the 2si/2 state by a horizontal green dotted line. 
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Fig. 5. Reduced supermomenta qK{f) (top panel) and full supermomenta 
Qnif) = Qnir) — n/r (bottom panel) for the 2si/2, «^ = 1 (red solid lines) and 
l'^3/2) = ~2 (blue dashed lines) pseudospin partner levels. 
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Fig. 6. Wave functions of the super symmetric partner states for the 2si/2 (red solid 
line) and 1^3/2 (blue dashed line) pseudospin partner levels. 
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